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STABILITY OF THE PARACHUTE AND HELICOPTER. 



By H. Batehan. 



INTRODUCTION. 



This paper was submitted to the National Advisory Committee for Aeronautics by 
Professor H. Bateman, of the California Institute of Technology, and its publication duly 
authorized by the committee as Technical Report No. 80. 

The mathematical theory of the stability of a parachute which is symmetrical with respect 
to a plane is very similar to the well-known theory of the stability of an airplane, but the values 
of the resistance coefficients are naturally different. When the parachute does not rotate, it 
may be compared with an airplane in a straight dive. There are oscillations corresponding to 
the longitudinal and lateral oscillations, and these are practically independent of one another. 
Thus, when the parachute swings like a pendulum as it descends, we have an oscillation corre- 
sponding to the pitching of an airplane, while there is another type of oscillation which corre- 
sponds to the Dutch roll. 

When the parachute is symmetrical about an axis, so that it has the form of a solid of rev- 
olution, there may be no banking in the oscillation corresponding to the Dutch roll, and the 
oscillation then becomes a simple swing like the one which corresponds to the pitching of an 
airplane. In this case a combination of two simple swings in perpendicular planes may give rise 
to a compound oscillation, which may be likened to the motion of a conical pendulum falling 
under gravity. In the general case the combination of the simple swing and Dutch roll gives 
rise to a compound motion of a complicated character. 



To discuss the matter analytically we shall write down the equations of motion in the 
notation used by G. H. Bryan 1 and S. Brodetsky.* We shall take the axis of y as axis of sym- 
metry and assume that in the steady state the parachute is falling with velocity V in a vertical 
direction, which coincides with the axis of y, and is rotating with angular velocity Q around the 
axis of y. This last assumption is made for the sake of generality, so as to cover the case of a 
rotating solid of revolution and to obtain a condition corresponding to some extent with the case 
of the helicopter where the longitudinal and lateral oscillations are not independent. The 
equations of motion may be written in the form 



ANALYTICAL DISCUSSION. 




+ (.Q + q)v)-r(V+v) =-Wsme-X, 
+p(V+v)-u(Q+q) = •- Wcos0sin0-Z, 
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"• Stability In Aviation. 

* The Toholcu Mathematical Journal, vol. 14, August, ISIS, p. 116. 
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Here Wis the weight of the parachute and. passenger, X, Y, Z, L, M, i\f the component 
air forces and couples referred to axes fixed in the parachute, B is the moment of inertia about 
the axis of symmetry, A the moment of inertia about a perpendicular axis through the center 
of gravity, both measured in gravitational units, U+u, v, w,_p, Q + q, r are component velocities 
and spins in the disturbed motion. To obtain the directions of the axes in space in a disturbed 
position of the parachute we suppose the parachute rotated first about the axis of y (from z to~x) 
through an angle then about the axis of z (from x to y) through an angle 6, and finally about 
the axis of x (from y to z) through an angle <j>. 

Let us write u + iw = a, p+ir=>h, $+£0=x where i denotes V - 1> then since u, w, p, r, 8, 
and $ are small in a small oscillation and the parachute is symmetrical about. the_ axis of y, wo 
may write to a first approximation 

d<j> d6 , dx 
P = Tt> r= Rt> i= dl 

X+iZ=>aE+bF 
L+iN=aJ+lK 

where the generalized resistance coefficients E, F, J, K are complex quantities. 
The above equations imply that if E= E t + iE t , etc., " 

Xj 1 a Z v =E 1] Z a = — Xy, = E 3) 

X P =Z T **F U Z v =—X T =F it 

iu=^w = J U iV a = -L^ =J„ 

Zp=iV r = Zi } JV,= -X r =K 2 , 

X v = 0, X<l=0, Z v =0, Z^O, 

It is on account of these relations that the analysis for the symmetrical parachute is 
simpler than that for an airplane in a straight drive. If, in fact, we neglect squares of small 

quantities and use D to denote the operatore ^ the first four equations may be written in the 
form 



Seeking a solution of the form a=a () e xt ' x=Xtfi u we obtain the period equation 

[-j(MG) +E ][f XS+ ( Z+ ^^) x ]" i; '[(-7^-*- F ') x+ F ] = 0 - This is a cubic nation 
with complex coefficients for the detennination of ,X. When Q = o it reduces to a cubic 
equation with real coefficients for we have the relations E t = o, F 1 = o, J 1 = o, jf£j = o. 

Remembering that in the general case^of a rigid airplane the oscillations about a state 
of steady motion are determined by an algebraic equation of the eighth degree, we infer that 
in the present degenerate case our real cubic-is a double factor of the octic and the periods and 
decrements are the same for simple swings in two perpendicular planes through the axis of 
symmetry of the parachute. 

When Q = o a root \=a+ip for which a is negative and jS positive implies the existence 
of a damped oscillation in which the center of gravity of the parachute describes a curve on a 
vertical cylinder in such a way that it moves in a counterclockwise direction when the axes 
are left handed. To see this we notice that u+_iw is of the form 

[cos 0OK,) +i sin 

where R and are real quantities. This equation indicates that ^ increases with t and since 
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the n-yis of 2 is to the left of the axis of x, the axis of y being 'downwards, this means that the 
center of gravity moves on the cylinder in the counterclockwise direction. 

The equation <t>+i8 = X 0 e ta+i M implies that at the same time the axis of symmetry pro- 
cesses around the vertical in the counterclockwise direction. If we disregard the vertical motion 
the character of the oscillation may be pictured by imagining a cone representing the para- 
chute to partly roll and partly slide on the outside of another cone whose axis is vertical. 

Let us call this a positive oscillation and use the term negative oscillation to denote one 
in which 0 is negative and the two motions take place in the clockwise direction. 

Either a positive or negative oscillation can be regarded as built up from two simple 
swings in perpendicular planes, the two swings having the same period and damping time, 
but differing in phase by a quarter period. By combining two such swings with an arbitrary 
difference in phase a type of elliptic oscillation is obtained. 

It should be noticed than an undamped positive oscillation may correspond to a new type 
of steady motion of the parachute. If such a motion exist it is characterized by v = o, q*=o 
and the fact that the cubic equation has two purely imaginary roots. 

Thus cubic equation may be written in the form 

X* +y\ 2 + xw\ + x = o 

*=V'^ = <X + TF> W "J~~W\ 1 TUT/' 
Let X 1 — be a factor of the cubic, then a and 9 are determined by the equations 

4a?-6+2o!y+xw=o 
—£ct8—yd+x=o. 

With the aid of these equations we can plot the curves t= constant in the (z, y) plane. 
Since 

Uci'+Say+xw) (y+2a)=x 
a curve £ = constant (a = constant) is a hyperbola. When a<=o it reduces to x = o, y = ^ and 
this is a boundary of the region of stability in the (z,y) plane. The condition for undamped 
oscillations is, moreover, either x = o or y=- % 

The asymptotes of the hyperbola t = constant are 

1 2a 
y+2a=~sad 2ay+xw + 4<^+—='0 

2a 

hence when 8= we must have x = oo y = oo . This means that we can not have a value of 8 

iff 

8a 

which is smaller than— — j consequently when the time of damping is given the period p can 

not be greater than a certain limiting value P. 

If a be given it is clear that P increases with w and so it is desirable that w should be 

made as large as possible. This can be done either by making V large or by making X T — ^^ T 

negative and fairly large. 

Now V is generally less than lg and ought not to be increased beyond this value, conse- 
quently the most hopeful method of improving the stability of the parachute is to make 

X T —^? T 1 negative. It is probable that in the Calthrop parachute this condition has been 

secured. 
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In Figure I the curves t => constant have been drawn for two different values of w. 
It will be seen that when w*=M the time of damping t is generally greater than the period. 
This is sometimes the case even when w = .6, but by choosing x and y in a suitable manner it is 
possible to make p greater than t. The importance of a fairly large value of w is thus manifest. 

In his study of the stability of the para- 
chute, Brodetsky treats the parachute as a 
circular disk attached to a weighted stick and 

finds from a condition equivalent to y > ^ that 

the center of gravity must lie between certain 
limits. . The center of gravity must naturally 
be fairly low, and there seems no way of 
avoiding this except perhaps by attaching a 
ring-shaped balloon to the rim of the disk. 
The upward thrust on the balloon and the 
downward force through the center of gravity 
would then have as their resultant a force 
acting through a point lower than the center 
of gravity, and this point would act as the 
center of gravity of the parachute in the 
usual theory. Unless it could be combined 
with a helicopter in some way, tho combi- 
nation of a parachute and toroidal balloon 
would be of theoretical interest only. 

Turning now to the case in which Q±0, 
we write 




/ 
<5 



F=iF„ + Q ft + ifj 



E=E 0 +Q (ci+ie,) 
J=iJ e + Q (j.+ih) 

The additional terras involving Q give 
the forces and couples arising from the Mag- 
nus effect or its inverse 1 and from related 



where 



4 6 6 

Fl<3. I. 

phenomena.* The cubic equation may now be written in the form 

, (iJ ^ + M±m.^(E 0 K a+ J.F.+ J a ^f)x + W. 

r= Ae 1 ±WJc 1 

9 ' """ ' ' 

WV . 

/= Efc + + J 0 f t + FJ, + — 

W WV B-A 

To study the effect of rotation upon the oscillations let o^ — i/3, be the roots 

of the equation $(X)=0, X 0 and being supposed negative and let <r 0 4~ir 0l <r t +ir u <r, + ir„ be 
the quantities which must.be added to X^^+iPu respectively, to give the correspond- 
ing roots of the equation ■&(/{.) =o. If o- and t are both positive it means that the effect of 
rotation is to diminish the period of a positive oscillation and to increase the time of damping f" 

1 For this effect see A. Laiay. Comptes RendusT Paris. ^oLlEJ&Ldf, pY^'voTIs! (1811)/^ !«!(.' - -- 
»For the effect of sideslip on a propeller see the Technical Reports ot the British Advisor j Committee lor Aeronautics 1912-13, 1MS-H, 
and 1918. 
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The effect on a negative oscillation is to increase both the period and the time of damping. 
The effect on a simple subsidence is to make it periodic and increase the time of damping. 
The effect of rotation when a and t are not both positive is easily inferred. 

In order that <s w <r 1; o- 3 may all be negative and the damping times for the three disturb- 
ances associated with i^, a 1 +ip 1 , o^— ifr all increased, it is necessary, but not sufficient, that 
QQ, QI, and Qj 2 should all be positive. Now N a = J t + Qj v hence it is necessary in the first 
place that the effect of the rotation should be to increase N^. The effect on the period is of 
some interest because in general the single period and time of damping associated with the 
two conjugate roots 0^+^ and e^—id^ will give rise to two different periods and two different 
times of damping; also a new positive or negative oscillation with a long period will arise from 
the simple subsidence associated with the root X a . The phenomenon of the division of one 
period into two is somewhat analogous to that which occurs in the Zeeman effect, especially 
as we may have a compound oscillation built up from positive and negative circular oscilla- 
tions of nearly equal period. The curves in Eigure II give the horizontal projection of the 
path of the center of gravity of the parachute or helicopter in two particular cases of com- 
pound oscillations of the above type. In the 
first case the two component circular oscilla- 
tions are undamped, in the second case they 
are both damped, the rate of subsidence being 
the same in each case. A complete discus- 
sion of the effect of rotation on stability is 
out of the question at present owing to our 
lack of knowledge of the values of the various 
resistance coefficients, but it may be worth 
while to ascertain the conditions for stability 
by finding the conditions that a cubic equa- 
tion with complex coefficients may have roots 
whose real parts are all negative. 

The conditions that a cubic equation with complex coefficients may have the real parts of all 
its roots negative. 

The method given by E. J. Eouth 1 for finding the conditions that an algebraic equation 
with real coefficients may have the real parts of all its roots negative may be extended to the 
case of an equation with complex coefficients. 

Let us consider the case of the cubic 

(x+aj (z+a,) (a + a s )^+ya? + g^+r=a?+(p l +ip 3 )^ + (2 1 +i2 ! )»+» , i+w a 

Let a^a^+ijh, ch = x i +iy 2 , a i = x s +iy s , and let the equation whose roots are — x 1 +iy lt 
-Xj+iy^ -Xi+iyt, be 

(x+p,)==x»+Px t +Qx+R. 

In Eouth's method the first step is to write iy for x and to separate the real and imaginary 
parts of the cubic, thus obtaining two expressions 

The process of finding the greatest common measure of f t (y) and f, (y) must now be carried 
out, the sign of the remainder being changed at each step just as in Sturm's theorem. In this 
way we obtain a series of polynomials whose first coefficients are 

i Loc.oit. 

J- 
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respectively! where 

H=p 1 p 3 &-q i 2 -p 1 r l +p 1 %, 

In order that-x 1} as,, and x, may all be positive.it is necessary and sufficient that p u H, 
and K should all be positive. The quantity K becomes zero when any one of the quantities 
x v Xt, x, becomes zero, it corresponds, therefore, to Eouth's discriminant in the case of the 
real quartic. In terms of x 1 , x,, Xg, y„ y v y„ the expressions for H and K are 1 

H = Ax t x, + Bx a x t + OxjXt, 
H'K^x&x, [B?- (W-nA) (mA-lB)-(mA-lB) (nB-mC)-(nB-mC) (JC-nA)] 

where 

A=(x t + X t ) @X l +X a + xJ + (ft-Vi)", 

B=(x i +x l ) (2x 2 +x,+x I ) + (y,-y 1 } , J 
O-ixt + x,) (2x i +x 1 +x i ) + (y 1 -y i ) 1 , 

It is not evident from these expressions that x u x„ and a;, are positive when p v R, and K 
are positive; 4 also B?K'\a of the eleventh degree in the quantities x u x,, x t , y v y v y t , and is not 
the simplest function of the coefficients which becomes .gejro when one-of the quantities z„ x st 
Xg is zero; consequently it will be worth while to find an alternative set of conditions which 
will make x,, and x a positive. 

Let us consider the_expressions 

T« («i+&> <«t+A> («i+A) C«i+ft) («»+A) (<*»+A) («.+&) («•+« («i+A) 

8 = [(«, + A) («• + A) + («. +A) («i + A) + («i + A) («a + A)] x 
[(«,+A) fa+A) + (<%+« K+A)lx 
[(«, +A) («. + A) + («. + A) («i + A) + («i + A) («» + A)l 

= 4[x,x,+x,x 1 +x 1 x J ] [{xj+Xg+ityg-y,)} {aJg+^+ityg-y,)} 
+ {x i +x i +t(y, -y t ) } {as, H-x, +i(y 1 -y 2 ) } + { & +t(y l -ft) } {*, +*, +i(y, -y,)}] X 
[{te+z,) {2 J +K,-»(y*-y»)}+{a5 J +a!.-t(yi-y s )} {^J-a^ + 

-S'=[(«j+A) fo+ft) + (<%+&) («i+A) + («i+A) C«.+A)]X 
[(<%+« («. + A3 + (««+ A) («i+A)+(«i+A) (o.+A)]x 
[(«.+A) (««+A) + (a.+A) («i+A) + («i+A) («»+A)] 
E^fo+Xg) + (yg-yj*] [4^+a!,) +Vh+x i y+ (y a -y,) 2 ]x 

[4x,(a\ + x,) + (an + x,) 1 + (y x - y,)*], 

Pl^S^ + Xj + X,. 

It will be noticed that the second and third factors of S are conjugate complex quantities 
and so their product is a positive quantity. If is clear then that if p„ S, and T are all pos- 
itive, the quantities 

Xj 4" Xj + Xg, XgXg + KfXi -f- XjSCfj XjXjXj 

are all positive.* Now let x v x^ and Xg be regarded as the distances of a point from the sides 
of an equilateral triangle. This is legitimate since Xi + Xj+Xg is positive. Keeping this 
last quantity constant a point can be used to represent any possible set of values of x^ x„ and 
Xg whose sum has this constant value. 

Now when XgXg +" x^ + x^ is positive the representative point lies within the circum- 
scribing circle of the triangle and when ay^x, is positive the representative point lies within 

> See note 1. 

i In the case of the quadrarlo Pi-Xi+Xt, 3-ZtXt l(zi+Xi) t +<l/i-9t)*l-PiPtfi-&+Pi'ti- Whan n and S are positive it Is evident that xi+tt 
and iii, are both positive and that consequently n and it are positive. 
• When these quantities are all positive S 1 Is also positive. 
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one of the four regions bounded by lines making an acute angle with one another. It is clear 
then that when ■ s rX i J r x t , xfa+xfa+xfa and xp#s t are all positive the representative point 
must lie within the triangle and x±, x^, Xg, be all positive. 

To express S and T in terms of the coefficients we notice that if two of the quantities 
ay, are interchanged, or if two of the quantities fi v |8 2 , j8 s are interchanged, S is trans- 

formed into S 1 . Hence we are justified in assu ming that there is an identity of type 

where 

A = 09, -ft) 

Comparing coefficients we find that £=—1. 
Writing* S in the form 

^=(2+e+?P-^ 1 -^2-^ 8 )C2+G+P?-^ 1 -^2-«A)(2+Q+?P-^i-a^ a -a 1 )S 8 ) 
= (.q+Q+pPHi+Q) + (2+Q+pP) (tfQ -t-PtqSqQ) - U 

where 

we assume that 

- U+^A=ArR+B(rPQ+RM) + °( rpa + R f)+DPQpq+Bp s P > + F(p t PQ + P i P q) 
Determining the unknown coefficients by putting 

(1) «1=0| Pi^o, 

(2) ft=o,/3,=o, 

(3) a,~afcft=o, 

(4) ft=&=/S„ 
we find that 

Hence finally 

S - ( 2 + Q) (SL+Q+pP)' + <a+Q+pP)&Q+P*q- 3 2 Q) -^tB +|(rPQ +Rpq) 
-TP*-Rtf-\pQpq-\&, 

A*-tff-4p>r + 18pqr-if-27i*) (P'Q* -4P*R +18PQR- 40 S -27P J ) 
Writing !F in the form 

r-W+PV + Qoi+B) C^+Paj' + Qaj+P) {aS+Paf + Qat+R) 
we easily find that 

T= r 8 + P s + Pr»2 + + W + 2^ - 2 QP^ - 2 2^ + 3Pr* + ZS?r + Rq* + r Q 3 - SRpqr - ZrPQR 
+ P i pr l + p*PR? + PV + p'P 3 + P^gr + pqPQR + P*Qqr + tfqQR + PQ'pr + pq>PR + F*Rq> 
+p 3 rQ 2 - 2P i Rpr - 2frPR - WPQ - 3R 2 pq - QRpr - grPR - 2q*QR - 2Q 2 gr. 

o 



A=~, B=\, C=-l, D~-\> E=o, F-o. 



